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We present a new method for treating the inner Cauchy boundary of a black hole spacetime by matching to a characteristic 
evolution. We discuss the advantages and disadvantages of such a scheme relative to Cauchy-only approaches. A prototype 
code, for the spherically symmetric collapse of a self-gravitating scalar field, shows that matching performs at least as well as 
other approaches to handling the inner boundary. 

Thus truncation of the interior spacetime at the ap- 
parent horizon does not affect the gravitational waves 
radiated to infinity. This is the physical rationale behind 
the apparent horizon boundary condition. 

There is a gauge ambiguity in the inner boundary de- 
fined by an apparent horizon which is associated with 
the choice of Cauchy foliation. Such an ambiguity is not 
associated with the event horizon. However, the event 
horizon is of no practical use in a Cauchy evolution since 
it can only be constructed in retrospect, after the global 
geometry of the spacetime has been determined. A bet- 
ter alternative is the trapping horizon Jl6{ , defined as 
the boundary of the spacetime region containing trapped 
surfaces. Here the reference to Cauchy hypersurfaces is 
dropped while retaining the quasilocal concept of trapped 
surfaces. Trapping horizons exist in any black hole space- 
time whereas the existence of apparent horizons is depen- 
dent on the choice of Cauchy foliation. 

In practice, the problem of locating trapped surfaces is 
partially solved in the process of setting initial data. For 
the 3-dimensional problem of two inspiraling black holes, 
there are several numerical approaches for determining 
appropriate initial Cauchy data ]l7t . An apparent hori- 
zon, when it exists, is a marginally trapped surface and 
lies on the trapping horizon. Once the initial Cauchy 
hypersurface cuts across a trapping horizon in this way, 
the scenario for pathological foliations is not present ini- 
tially; and a reasonable choice of lapse should guaran- 
tee that future Cauchy hypersurfaces continue to contain 
that component of the apparent horizon. However, in 
the two black hole problem, besides the two disjoint ap- 
parent horizons present initially, an outer apparent hori- 
zon (surrounding them) is expected to form at a later 
time. Finding and locating this outer apparent horizon 
can make the computational problem enormously eas- 
ier by using it as the new inner boundary at this stage. 
Excellent progress has been made in designing apparent 
horizon finders and trackers for this purpose. However, it 
is not known what lapse condition on a Cauchy foliation 
would guarantee that an outer apparent horizon form at 
the earliest possible time. 

Besides these geometrical issues there are a number 
of serious computational difficulties in implementing an 
apparent horizon boundary condition. In order to ob- 
tain gravitational waveforms, the computational domain 



I. INTRODUCTION 

In many physical systems, boundary conditions are 
both the most important and the most difficult part of 
a theoretical treatment. In computational approaches, 
boundaries pose further difficulties. Even with an an- 
alytic form of the correct physical boundary condition 
in hand, there are usually many more unstable numer- 
ical implementations than stable ones. Nowhere is the 
boundary problem more acute than in the computation 
of gravitational radiation produced in the coalescence of 
two black holes. In order to avoid the topological com- 
plications introduced by the black holes, the proposed 
strategy for attacking this problem, initially suggested by 
W. Unruh B, is to excise an interior region surrounded 
by an apparent horizon. These are uncharted waters and 
there are many different tactics that can be pursued to 
attain an apparent horizon boundary condition. P-|l3|] 
One common feature of all current approaches to this 
problem is the use of a Cauchy evolution algorithm in 
the interior region bordering the apparent horizon. In 
this paper we present an alternative tactic based upon 
a characteristic evolution in that inner region; and we 
present a simple model of its global implementation. 

In order to provide orientation, we begin with a syn- 
opsis of the apparent horizon boundary condition and its 
computational difficulties. An apparent horizon is the 
boundary of the region on a Cauchy hypersurface con- 
taining trapped surfaces ]l4|j . This explicit reference to 
a Cauchy hypersurface in the definition gives an appar- 
ent horizon an elusive nature. Indeed, there are Cauchy 
hypersurfaces in the extended Schwarzschild spacetime 
which come arbitrarily close to the final singularity but 
do not contain an apparent horizon |l5[ . There is strong 
reason to believe that the same is true in any spheri- 
cally symmetric black hole spacetime. On the other hand, 
when they exist, apparent horizons are useful spacetime 
markers because they must lie inside the true event hori- 
zon |14j . Consequently, signals cannot propagate causally 
from the apparent horizon to future null infinity X + . 
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must cover a time interval of the order of several hundred 
M in the exterior region whereas typically a singularity 
forms on a time of order M in the region close to the 
apparent horizon. Thus a slicing which avoids the sin- 
gularity for several hundred M will necessarily develop 
coordinate singularities. In addition, the inner boundary 
traced out by an apparent horizon is generically space- 
like (at best lightlike). Thus if the coordinates defin- 
ing the numerical grid were to remain constant in time 
on the boundary ( "apparent horizon locking" ) then the 
coordinate trajectories would have to be superluminal. 
While horizon locking works in the spherically symmet- 
ric case it is difficult to implement in a Cartesian 
3-dimensional grid. The alternative is to let the apparent 
horizon move through the coordinate grid. At the same 
time, the location of the apparent horizon must be de- 
termined by solving an elliptic equation or an equivalent 
extremum problem. The requirements on the grid are 
further complicated when the black hole is spinning. On 
top of all these difficulties, the computational techniques 
must ensure that the strong fields inside the apparent 
horizon boundary do not severely leak into the exterior 
due to finite difference approximations. Causal differenc- 
ing Q and algorithms based upon a strictly hyperbolic 
version of the initial value problem jL8| have been pro- 
posed to avoid this. However, no 3-dimensional Cauchy 
code has yet been successful in evolving a Schwarzschild 
black hole. 




FIG. 1. A schematic of how matching to ingoing null cones 
could be used with two black holes. The inner Cauchy evolu- 
tion is matched at an outer world tube C to a null evolution 
on outgoing null cones; and at two interior world tubes A and 
B, to null evolutions on ingoing null cones. The evolutions on 
the ingoing null cones stop at the apparent horizons (dotted 
lines) which surround the two black holes. 

It is clear that the 3-dimensional coalescence of black 
holes challenges the limits of computational know-how. 
We wish to present here a new approach for excising an 
interior trapped region which might provide enhanced 
flexibility in tackling this important problem. In this ap- 
proach, we locate the interior boundary of the Cauchy 
evolution outside the apparent horizon. Across this in- 
ner Cauchy boundary we match to a characteristic evolu- 
tion based upon an ingoing family of null hypersurfaces. 



It is the inner boundary condition for the characteris- 
tic evolution which is then given by a null hypersurface 
version of the apparent horizon boundary condition. In 
the case of two black holes, the inner boundary would 
consist of two disjoint topological spheres, chosen so that 
their inner directed null normals are converging. FIG 
[l] provides a schematic picture of the global strategy. 
Two disjoint characteristic evolutions, based upon ingo- 
ing null hypersurfaces, are matched across worldtubes A 
and B to a Cauchy evolution of the shaded region. The 
interior boundary of each of these characteristic evolu- 
tions borders a region containing trapped surfaces. The 
outer boundary of the Cauchy region is another world- 
tube C, which matches to an exterior characteristic evo- 
lution based upon outgoing null hypersurfaces extending 
to null infinity. 

This strategy offers several advantages in addition to 
the possibility of restricting the Cauchy evolution to 
the region outside the black holes. Although, finding a 
marginally trapped surface on the ingoing null hypersur- 
faces remains an elliptic problem, there is a natural ra- 
dial coordinate system (r, 9, <fi) to facilitate its solution. 
However it is also possible to locate a trapped surface 
on the ingoing null hypersurface by a purely algebraic 
condition. Since this trapped surface (when it exists) 
lies in the region invisible to X + it can be used to re- 
place the trapping horizon as the inner boundary. In 
either case, moving the black hole through the grid re- 
duces to a 1-dimensional radial motion, leaving the an- 
gular grid intact and thus reducing the complexity of 
the computational masks which excise the inner region. 
(The angular coordinates can even rotate relative to the 
Cauchy coordinates in order to accommodate spinning 
black holes). The chief problem of this approach is that 
a caustic may be encountered on the ingoing null hyper- 
surface before entering the trapped region. This is again 
a problem whose solution lies in choosing the right initial 
data and also the right geometric shape of the two-surface 
across which the Cauchy and characteristic evolutions are 
matched. There is a great deal of flexibility here because 
of the important feature that initial data can be posed 
on a null hypersurface without constraints. 

The strategy of matching an interior Cauchy evolu- 
tion to an exterior outgoing characteristic evolution has 
been described |l9|-pl|] and implemented to provide a 
computational Cauchy outer boundary condition in var- 
ious cases, ranging from 1 and 2 dimensional simula- 
tion s |22| |25|] to 3-dimensional simulations that include 
T + |26| , |27[ . A slight modification allows changing an out- 
going null formalism (and its evolution code) to an ingo- 
ing one. This is briefly reviewed in Sec. |J. By matching 
Cauchy and characteristic algorithms at both an inner 
and outer boundary, the ability to include I + facilitates 
locating the true event horizon while excising an interior 
trapped region. In Sec. Ill, we discuss the problem of 



locating trapped surfaces on an ingoing null hypersur- 
face. In Sec. IV, we present an implementation of these 



ideas to the global evolution of spherically symmetric, 
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self gravitating scalar waves propagating in a black hole 
spacetime. In this case, the performance of the matching 
approach equals that of previous Cauclw-only schemes 
that have been applied to this problem p-|5|,p8| . 



consistent because contains a free integration constant 
which can be chosen to be complex (as long as it leads 
to a real metric). In order to see how this works con- 
sider the outgoing version of the null hypersurface equa- 
tions fl3lJ33|: 



II. CAUCHY-CHARACTERISTIC MATCHING 

A. The null formalism 

We introduce a unified formalism for coordinates based 
upon either ingoing or outgoing null hypersurfaces. Let 
w label these hypersurfaces, x A (A = 2,3), be labels for 
the null rays and r be a surface area distance. In the 
resulting x a — (w, r, x A ) coordinates, the metric has the 
Bondi-Sachs form p^ , pp[ 

ds 2 — g W wdw 2 + 2g wr dwdr + 2g w Adwdx A 

+g AB dx A dx B , (2.1) 

where det(gA B ) = r 2 det(qA B ) — r 2 q, with qAB a unit 
sphere metric. In the outgoing case, writing w = u, it is 
convenient to express the metric variables in the form 

ds 2 = - (e W ^ - r 2 h AB U A U B ^j du 2 - 2e 2f3 dudr 

- 2r 2 h AB U B dudx A + r 2 h AB dx A dx B , (2.2) 

where h AB h B c = $c- This yields the standard outgo- 
ing null coordinate version of the Minkowski metric by 
setting = U A = V — r = Hab — Q_ab = 0. In the 
ingoing case, writing w = v, the only component of the 
Minkowski metric which differs is g vr = —g U r- This can 
be effected by the substitution 
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0.r = T^ rhAChBDh AB,rhcD,r 

16 

+2nrT rr , (2.5) 
(r*e- 2 Ph AB U B ), r = 2r 4 (r- 2 /^) r - r 2 h BC D c h AB , r 

+16nr 2 T rA (2.6) 
2e- 2l3 V r = 11- 2D A D A - 2D A 0D A 
+r- 2 e-^D A (r 4 U A ) tr 

-\r*e-^h AB U A U B 

+87rr 2 (T-g AB T AB ), (2.7) 

where Da is the covariant derivative and 1Z th e cu rvature 
scalar of the 2-metric Iiab- The /3-equation (2.5) allows 
the substitution ( |2.3|) to be regarded as a change in inte- 
gration co nsta nt . Th en carrying out this substitution in 
equations (2.5)-(2.7) leads to the ingoing version of the 
null hypersurface equations: 



13 -> + m/2 



(2.3) 



in the outgoing form of the metric. 

The substitution (2.3) can also be used in the curved 
space case to switch from outgoing to ingoing coordi- 
nates, in which case it is equivalent to an imaginary shift 
in the integration constant for the Einstein equation de- 
termining (see equation (2.5) below). This leads to the 



ingoing version of the metric 



,V 



ds 2 = e 213 - + r 2 h AB U A U B dv 2 + 2e 2p dvdr 



2r J, h AB U B dvdx A + r 2 h AB dx A dx 



(2.4) 



Of course, at a given space-time point the values of the 
coordinates r and x A and the metric quantities 0, U A , 
V and Ka B are not the same in the ingoing and outgoing 
cases; but since we do not consider transformations be- 
tween outgoing and ingoing coordinates there is no need 
to introduce any special notation to distinguish between 
them. 

This same substitution also provides a simple switch 
from the outgoing to the ingoing version of Einstein equa- 
tions G a /3 = 87rT Q/ 3 written in null coordinates. This is 



0,r = 



^rh AC h BD h A B.rhcD,r 



-(r 4 e- 2 ^h AB U B ), r 
-2e- 2 ?V r 



1 

f6 

+2irrT rr , 
2r 4 (r- 2 0,A), r ~ 

+16nr 2 T r A 
K 

■ 20 D A (r 4 U A 

-4/3 



r 2 h BC D c h AB , r 



(2.8) 

\ B ,r 

(2.9) 



2D A D A - 2D A 0D A 



-r 2 e 

1 4 - 

--re 



IiabU a U b 



-Mr z {T-g AB T AB ). 



(2.10) 



This formal substitution also applies to the dynamical 
equations and provides a simple means to switch between 
evolution algorithms based upon ingoing and outgoing 
null cones. 

As we have already noted, although the same coor- 
dinate labels r and x A are used for notational simplic- 
ity in both the outgoing metric ( |2.2| ) and the ingoing 
metric (2.4), they represent different fields. An excep- 
tion occurs for spherical symmetry where the the surface 
area coordinate r can be defined uniquely in terms of the 
same 2-spheres of symmetry used in both the ingoing and 
the outgoing coordinates. In this case, the spacelike or 
timelike character of the r = const hypersurfaces is con- 
sistent under the substit utio n (|2.3|) be cause the change 
involved in going from ( |2.7| ) to (| 2 . 1 C ) implies that V 



changes sign in switching from outgoing to ingoing co- 
ordinates. As a result we obtain a consistent value for 
g a Pr tU r^ = ±e~ 2/3 V/r, with the + (— ) sign holding for 
outgoing (ingoing) coordinates. 
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In the absence of spherical symmetry, the surface area 
coordinate r used in the Bondi-Sachs formalism has a 
gauge ambiguity associated with the changes in ray la- 



bels 



y (x ), under which it transforms as a scalar 



density. On any null hypcrsurface with a preferred com- 
pact spacelike slice So, this coordinate freedom in r may 
be fixed by requiring that r = const on So- This then 
determines a unique r = const foliation on either the in- 
going or outgoing null hypersurface emanating from Sq. 



B. Matching 

Cauchy-characteristic matching can be used to replace 
artificial boundary conditions which are otherwise neces- 
sary at the outer boundary of a finite Cauchy domain. 
The exterior characteristic evolution can then be ex- 
tended to null infinity to form a globally well-posed initial 
value problem. In tests of nonlinear 3-dimcnsional scalar 
waves, Cauchy-characteristic matching dramatically out- 
performs the best available artificial boundary condition 
both in accuracy and computational efficiency ]26| , p7| ]. 

We now describe how this matching strategy can be 
used at the inner boundary of a Cauchy evolution which 
is joined to an ingoing null evolution. On the initial 
Cauchy hypersurface, denoted by time to, let So be a 
(topological) 2-sphere forming the inner boundary of the 
region being evolved by Cauchy evolution. Let W rep- 
resent the future evolution of So under the flow of the 
vector field t a — an" + /? a , where n a is the unit vector 
field normal to the Cauchy hypersurfaces t — const and a 
and f3 a are the lapse and shift. Given the initial Cauchy 
data on to, boundary data must be given on the world 
tube W in order to determine its future evolution. The 
Cauchy hypersurfaces foliate this world tube into spheres 
St- The boundary data is obtained by matching across 
W to an interior null evolution based upon the null hy- 
persurfaces v = const emanating inward from from the 
spheres St- The evolutions are synchronized by setting 
v = t on W. 



FIG. 2. The shaded region is the domain of dependence of 
initial data which is given on a Cauchy slice io and an ingoing 
null slice vq. Q is the apparent horizon. W is the matching 
world tube. 

The combination of initial null data on vq and initial 
Cauchy data on to determine the future evolution in their 
combined domain of dependence as illustrated in FIG [|. 
In order to avoid dealing with caustics we terminate the 



null hypersurfaces vt on an inner boundary Q whose loca- 
tion is determined by a trapping condition, as discussed 
below. This inner boundary plays a role analogous to an 
apparent horizon inner boundary in a pure Cauchy evo- 
lution. The region inside Q is causally disjoint from the 
domain of dependence which is evolved from the initial 
data. In the implementation of this strategy in the model 
spherically symmetric problem of Sec. IV, we describe in 
detail how data is passed back and forth across W to sup- 
ply an outer boundary value for the null evolution and 
an inner boundary value for the Cauchy evolution. In 
the remainder of this section, we discuss how some of the 
key underlying issues might be handled in the absence of 
symmetry. 

In order to ensure that an inner trapping boundary 
exists it is necessary to choose initial data which guar- 
antees black hole formation. Such data can be obtained 
from initial Cauchy data on to for a black hole. However, 
rather than extending the Cauchy hypersurface inward 
to the apparent horizon it could instead be truncated at 
an initial matching surface So located sufficiently far out- 
side the apparent horizon to avoid computational prob- 
lems with the Cauchy evolution. The initial Cauchy data 
would then be extended into the interior of So as null data 
on vo until a trapping boundary is reached. Two ingredi- 
ents are essential in order to arrange this. First, So must 
be chosen to be convex, in the sense that its outward null 
normals uniformly diverge and its inner null normals uni- 
formly converge. Given any physically reasonable mat- 
ter source, the focusing theorem then guarantees that the 
null hypersurface vq emanating inward from So continues 
to converge until reaching a caustic. Second, initial null 
data must be found which leads to trapped surfaces on 
Vq before such a caustic is encountered. The existence 
of trapped surfaces depends upon the divergence of the 
outward normals to slices of vo- Given the appropriate 
choice of So the existence of such null data is guaran- 
teed by the evolution of the extended Cauchy problem. 
However, it is not necessary to actually carry out such a 
Cauchy evolution to determine this null data. It is the 
data on So which is most critical in determining whether 
a black hole can form. This can be phrased in terms 
of the trapping gravity of So, introduced below. In the 
spherically symmetric Einstein-Klcin-Gordon model (see 
Sec. IV), if So has sufficient trapping gravity to form 
a trapped surface on vo in the absence of scalar waves 
crossing vq then a trapped surface also forms in the pres- 
ence of scalar waves. In the absence of symmetry, this 
suggests that given appropriate initial Cauchy data for 
horizon formation that the simplest and perhaps physi- 
cally most relevant initial null data for a black hole would 
correspond to no gravitational waves crossing vo- The 
initial null data on vo can be posed freely, i.e. it is not 
subject to any elliptic or algebraic constraints other than 
continuity requirements with the Cauchy data at iSo- In 
the vacuum case, the absence of gravitational waves in 
the null data has a natural (although approximate) for- 
mulation in terms of setting the ingoing null component 
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of the Weyl tensor to zero on Vq. 

Key to the success of this approach is the proper trap- 
ping behavior, i.e. the convergence of both sets of null 
vectors normal to a set of slices of «o located between 
the caustics and the matching boundary. By construc- 
tion, the ingoing null hypersurface Ao, given by v = Vq, 
is converging along all rays x A leaving the initial slice 
So coordinatized by r = Rq. (Here (v,r, x A ) are ingoing 
null coordinates.) In order to investigate the trapping of 
Ao we must determine the divergence of slices S of Ao 
defined by r — R{x A ). Let n a be tangent to the gen- 
erators of Ao, with normalization n a — —g V rV,a- Then 
n Q w jQ = n ^^, = and n Q r Q = —1. Let l a be the 
outgoing normal to 5, normalized by n a l a = — 1. Then 



1 



-lV, a + r, a 



R,AX A a 



where 



(2.11) 



III. TRAPPING 



A slice of Ao is trapped if @i < 0. I n te rms of the 
ingoing Bondi metric variables defined in (2.4), equation 
( [2.16 ) gives 



r 2 e 2/3 



e, 



1 



V-^[Vq(e 2l3 h AB R B -r 2 U A )U 



~ r{r^e 2l3 h AB )^R A R^ + r 2 R, A U A . (3.1) 



Setting Qi = in equation ( p.l[ ) gives a 2-dimensional 
Laplace equation for the function R(x A ) which locates a 
marginally trapped surface A. Such a surface lies on a 
trapping horizon and is (a component of) the apparent 
horizon of any Cauchy hypersurface which contains it. 

For a marginal surface to lie on an outer trapping hori- 
zon its trapping gravity, defined as p6[ 



-g vr (g rr ~ 2g rA R iA + g AB R. A R, B ). (2.12) 



Its contravariant components are 



r 



i r = \g rr - \g AB R,AK 



and 



g R.b- 



(2.13) 



(2.14) 



(2.15) 



Let lp — g'p + n a l/3 + l a np be the projection tensor 
into the tangent space of S and define 7 Q/3 and 7 Q/ 3 by 
raising and lowering indices with g a p. Its contravariant 
components are j av = 0, j AB = g AB , j rA = g AB R,B 
and 7 rr = g AB R, a R : b', and its covariant components 
are 7 Qr = 0, 7,45 = g AB , -y v A = g vr R,A + g v A and ^ vv = 

g AB {gv r R,A + gvA)(g vr R,B + gvs)- 

The outward divergence of S is given by 9; = 
2^ Q / 3 V Q //3. (The conventions are chosen so that = 
2/r for a r — const slice of an outgoing null cone in 
Minkowski space.) Then a straightforward calculation 
yields 

1 1 a vr 

-0, = -g rr - ^[^g AB (g V rR,B+g V B)U 

AB\ 



- — {rg vr g AB ), r R,AR 
r 

-g vr R, B (9 AB 9vA), r , 



4c 



(3.2) 



(2.16) 



must be real and positive, so that k > 0. The trapping 
gravity generalizes the concept of the surface gravity of 
an event horizon to trapping horizons. In our coordinate 
system, n a d a — ~d r . Thus if the surface r — R(x A ) is 
marginally trapped then positive trapping gravity implies 
that the surface r — R (x A ) — Ar is trapped for small Ar. 
We use equation ( |3.2|) to generalize the definition of trap- 
ping gravity to an arbitrary slice of a inwardly converging 
null hypersurface. Then slices of positive trapping grav- 
ity tend toward trapping as r decreases. However, in 
general, there seems to be no purely local criterion which 
guarantees the existence of a trapped surface before en- 
countering a caustic as r — > 0. 

In the special case of a spherically symmetric slice of a 
spherically symmetric null cone, the Laplace equation for 
a marginally trapped slice reduces to the algebraic con- 
dition that V = 0, which is satisfied where the r = const 
hypersurface becomes null. The vacuum Schwarzschild 
metric in ingoing null coordinates (which are equivalent 
to ingoing Eddington-Finkelstein (IEF) coordinates) is 
given by j3 = 0, V = -(r-2M), U A = and h AB = q AB - 
In this case, V = determines the location of the event 
horizon (which coincides with the apparent horizon) and 
the surface gravity reduces to k = (4Af) _1 . In the non- 
vacuum spherically symmetric case, V = determines 
the apparent horizon. 

In the absence of spherical symmetry, 0; vanishes on 
a slice of the form R = const at points for which Q = 0, 
where 



which is to be evaluated on S after the derivatives are 
taken. The divergence of the generators tangent to Ao is 
given by 6„ = j a0 V a n p . Then 6„ = -2/r so that A" 
is converging in accord with our construction. 



Q 



-V- 



v 7 ? 



(Vtu A ), A . 



(3.3) 



We will refer to the largest r = const slice of Ao on 
which Q < as a "Q-boundary" Q, relative to S . (S 
enters here because it provides the reference for r = const 
slices). Such a slice is everywhere trapped or marginally 
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trapped so that the Q-boundary provides a simple alge- 
braic procedure for locating an inner boundary inside an 
event horizon. A Q-boundary, when it exists, will always 
lie inside (smaller r) or tangent to the trapping hori- 
zon. Let r = Rq describe the Q-boundary and let R m in 
(Rmax) be the smallest (largest) value of r on the trap- 
ping horizon. Then on the trappi ng ho rizon R : a = at 
R m in and R m ax so that equation ( 2.16 ) implies Q > at 
Rmin (and Q < at R max )- Consequently, Rq < R min , 
with equality holding only when D 2 R = at R m in- 

There are thus two possible strategies for positioning 
an inner boundary, both of which ensure that the ig- 
nored portion of spacetime cannot causally effect the ex- 
terior spacetimc: (I) Use the trapp ing horizon, in which 
case the 2D elliptic equation ( 2.16 ) must be solved on a 
sphere in order to determine its location; or (II) use the 
Q-boundary which is determined by a simple algebraic 
condition. 

Strategy (I) is similar to approaches used to locate an 
apparent horizon on a Cauchy hypersurface. The ad- 
vantage in the null cone case is that there is a natural 
radial coordinate defined by the coordinate system to re- 
duce the elliptic problem to 2 angular dimensions and 
to define a mask for moving the excised region through 
the computational grid. Strategy (II) carries no essential 
computational burden since the quantities V and U are 
obtained by means of an inward radial integral as part of 
the evolution scheme. One merely stops the integration 
when the inequality defining the Q-boundary is satisfied. 
Thus strategy (II) is preferable unless either caustics or 
singularities appear before reaching the Q-boundary. It 
is easy to choose black hole initial data so that the Q- 
boundary and trapping horizon agree at the beginning of 
the evolution but whether the Q-boundary will move too 
far inward to be useful is a critical question which would 
depend upon choices of lapse, shift and the geometry of 
the matching world tube. Further research is necessary 
to decide if strategy (II) is viable on geometric grounds 
in a highly asymmetric spacetime. 



Our purpose here is to present a spherically symmetric 
model demonstrating the feasibility of a stable global al- 
gorithm based upon three regions which cover the space- 
time exterior to a single black hole. FIG || illustrates 
1-dimensional radial geometry. The innermost region 
is evolved using an ingoing null algorithm whose inner 
boundary Q lies at the apparent horizon and whose outer 
boundary Rq lies outside the black hole at the inner 
boundary of a region evolved by a Cauchy algorithm. 
Data is passed between these regions using a matching 
procedure which is detailed below. The outer bound- 
ary i?i of the Cauchy region is handled by matching to 
an outgoing null evolution. The details of the outgoing 
null algorithm |33j | and of the Cauchy evolution || are 
not discussed since they have been presented elsewhere. 
We will discuss the matching conditions since they differ 
from those used previously EfJ due to different choices of 
gauge conditions. We will also present the field equations 
since they are important for understanding the matching 
procedure. 

The Cauchy evolution is carried out in IEF coordi- 
nates. The metric in this coordinate system is 

ds 2 = a 2 (20 - l) dP + 2a 2 0dtdr + d 2 dr 2 + r 2 dfl 2 . 

(4.1) 

The set of equations used in the evolution are 



_4gn 



-a 2 



1-0) K 



dt3 



(4.2) 



(4.3) 



(4.4) 



IV. COLLAPSE OF A SPHERICALLY 
SYMMETRIC SCALAR WAVE 



+ 



1-0 



K% = 0K 6l e +a(l-0)K e e {K\ + 2K 



0' 



(4.5) 





FIG. 3. This figure shows the two matching world tubes 
and the three coordinate patches. 



$ = 



(/3$+ (1-0) n)' : 



f z (0ii + 1 - $ 



(4.6) 



(4.7) 



where the over-dot represents partial with respect to t, 
prime denotes partial with respect to f and the scalar 
field variables are defined by 



n 



1 



1-0 



(4.8) 
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The outgoing-null metric is 



e 2l3 —du 2 - 2e 2l3 dudr + r 2 dfl 2 , 



(4.9) 



The outgoing hypersurface equations (2.5) - (2.7) reduce 
to 



/? r = 2nrq 
— e" 



V r - ^ 



(4.10) 
(4.11) 



with U A = and h AB = qAB, and the outgoing version 
of the scalar wave equation is 



2M, r „ = -{rV<t>, 
r 



(4.12) 



An ingoing null evolution algorithm can be obtained 
from the outgoing algorithm by the procedure described 
in Sec. [n| Given the ingoing null metric 



ds 2 = e 2p —dv 2 + 2e 2f3 dvdr + r 2 dU 2 



(4.13) 



an independent set of equations are the hypersurface 
equations 



T = 2-nrq 



=2/3 



and the scalar wave equation 

2(n/>), r „ = -{rV<f>, r ),r 
r 



(4.14) 
(4.15) 



(4.16) 



Given data for (j> on Ao and on the world tube IZo, defined 
by r — i?oi an d integration constants for (3 and V on 1Z 
evolution proceeds to null hypersurfaces J\f v , defined by 
v = const, by an inward radial march along the null rays 
emanating inward from IZq. 



A. Initial data 

The initial data consists of a Schwarzschild black hole 
of mass M which is well separated from a Gaussian pulse 
of (mostly) ingoing scalar radiation. Initially there is 
no scalar field present on the ingoing-null patch, so the 
initial data there is simply 



= 0, 
[3 = 0, 

V = 2M - r. 



(4.17) 
(4.18) 
(4.19) 



Similarly, initially there is no scalar field on the 
outgoing-null patch, so we have = 0. The values for (3 
and V are determined by matching to the Cauchy data 
at the world tube (see below) and integrating the hyper- 
surface equations (4.14) and (4.15). 

In the Cauchy region, the scalar field is given by 



Arexp 



n 



(f — c) c I a 



1 d(r-c) 
r o~ d 



2-/3 d(f-c) 



d-l 



r (1-/3 



(4.20) 
(4.21) 

(4.22) 



where A, c, d, and a are scalars representing the pulse's 
amplitude, center, shape, and width, respectively. The 
geometric variables are initialized using an iterative pro- 
cedure, as detailed in |^). 



B. Matching conditions 

Since the IEF coordinate system is based on ingoing 
null cones, it is possible to construct a simple coordinate 
transformation which maps the IEF Cauchy metric to 
the ingoing null metric, namely 



t = v — r, 



(4.23) 



This results in the following transformations between the 
two metrics. 



= 1/2 log 
2/3-1 



M) 



V = r- 



1 

V - 



-/3 

■ r 



13 V + 2r 
a = ePy/V/r + 2. 



(4.24) 

(4.25) 

(4.26) 
(4.27) 



The extrinsic curvature components can be found using 
only the Cauchy metric, 



K\ = 



K r r = 



M) 



5 1-/3 



a a 



(4.28) 
(4.29) 



Note that these transformation equations are valid ev- 
erywhere in the spacetime, not just at the world tube. 

The matching conditions at the outer world tube are 
more complicated. Both the Cauchy and characteristic 
systems share the same surface area coordinate r but 
there is no universal transformation between their corre- 
sponding time coordinates. However, we can construct 
a coordinate transformation which is valid everywhere 
on the world tube. To do this, we start with a general, 
differential coordinate transformation, whose unknown 
function F is to be determined from the matching con- 
ditions: 



du = F(t)dt — df, dr = dr. 



(4.30) 
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To keep the notation simpler, we will write the Cauchy 
metric as 



ds = goodt + 2goxdtdf + g\\df 



- 2 ' ~ 2 dQ 2 



(4.31) 



Inverting (4.3C) and substituting, we get 



ds 2 = g m -^du 2 + 2 (J^ga + J^doo^j dudr 

+ (511 + JtfOi + ^5oo) dr 2 + rW. (4.32) 



On the world tube, we require u = t, sothat du = dt. 
Thus, we set du — du/F. For the metric, we get 

ds 2 = goodu 2 + 2 (g 01 + -^500^ dudr 

3ii + j9oi + ^.9oo) dr 2 + r 2 dn 2 . (4.33) 



The condition that the r-direction be null implies that 

(4.34) 



2 1 
.911 + p9oi + ^2.9oo = 0. 



Upon substitution of the IEF metric functions, this de- 
termines that 



F = 1-2(3. 
The matching conditions are then 



= 1/2 log 

V = r =- 

V-r 



(1-,) 



V - 2r 



a = e l3 ^2-V/r. 



(4.35) 

(4.36) 

(4.37) 

(4.38) 
(4.39) 



C. Finite difference implementation 

As is typical with finite difference calculations, the 
continuum functions are discretized spatially and placed 
on grids with TV points. In our case we have three re- 
gions. The inner-null variables are placed on grids with 
Ni points, the Cauchy variables on grids with N c points, 
and the outer-null variables on grids with N a points. On 
each grid, the spatial index i runs from 1 to TVj, N c or 
N Ql respectively, with 1 representing the smallest r value. 




3 

□ 

□ 
□ 
□ 



n+2 
n+1 
n 

n-1 



FIG. 4. This diagram shows the finite difference grids 
around the inner world tube. The squares represent the 
Cauchy grid, while the circles are the null grid. The diag- 
onal lines are the ingoing null cones, and the vertical line Ro 
is the inner world tube. The innermost points of the Cauchy 
grid lie on the world tube, while the null grid extends outside 
the world tube. Notice that the grids align in space and time. 



In addition to the spatial discretization, each function 
needs two or more time levels. While both the Cauchy 
and null evolution schemes use only two time levels, we 
keep an extra level in each to facilitate the matching. 
FIG ^ shows how the finite difference grids match at the 
inner world tube. The two grids are aligned in both r and 
t. This means no interpolations are necessary. The world 
tube is at i — 1 on the Cauchy grid and i = Ni — 1 on 
the ingoing null grid. The Cauchy variables need bound- 
ary values on time level n + 1 at i = 1. The metric 
values come from the null variables at level n+1, with 
i = Ni — 1, using the transformation equations (4.2£) 



and (4.27). These relations are algebraic and straight- 
forward to implement. Boundary values for th e ext rinsic 
curvature components come from equations (4.28) and 
(4.29). K e g is computed algebraically and K r r is com- 



puted using second-order, centered-in-time, forward-in- 
space derivatives in the Cauchy grid. 

The transformation of the scalar field requires trans- 
forming derivatives between the two coordinate patches. 
At the world tube, the relationships among the deriva- 
tives are 

d r = df- d t (4.40) 
dr = d r + d v (4.41) 
d v = d t . (4.42) 

These lead to the following equations for $ and II: 

9%- 



1 n n+2 
$ n+l _ J S \ 1 



n+1 
HNi-2 



2dv 



2dr 



n+1 
«2 



(4.43) 
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and 



n n+i 



9'k- 



n+l 



2dv 



n+l 



(4.44) 



The null variables need boundary values at n + 2, i = 
Ni. The metric values come from the Cauchy variables at 
level n, with i = 2, using Q4.24D and ( |4.25| ). The evolution 
equation for g at the world tube is 



9t 2 



at 2 



+ 



(4.45) 



n+2 
n+l 
n 

n-1 



N -2 

c 

□ 
□ 
□ 
□ 




FIG. 5. This diagram shows the finite difference grids 
around the outer world tube. The squares represent the 
Cauchy grid, while the circles are the null grid. The diagonal 
curves are the outgoing null cones, and the vertical line Ri 
is the outer world tube. The outermost points of the Cauchy 
grid he on the world tube, while the null grid extends inside 
the world tube. Notice that the the grids align in time only 
on the world tube, but in space at the world tube and just 
inside. 



The situation at the outer world tube is shown in FIG 
|| Here, the grids align in space at two values of r but in 
time only at the world tube. The Cauchy metric bound- 
ary values at i = N c come directly from the nul l vari ables 
at i = 2 using the transformation equations (4.38) and 
Since the grids don't align in time as they do 



Q4.39D 



at the inner world tube, we use a different procedure for 
the scalar field boundary values. We use the derivative 
transformation 



d r = 



1 

1 - 2/3' 



(4.46) 



the world tube. 

and 



to get an evolution equation for <p at 
We then set $ and II using their definitions (4 
backward second-order derivatives. 

The boundary values for the null variables must be 
interpolated in time using the Cauchy variables at n and 
n—1. Given values for a function / at time levels n and 
n— 1, we can get its value f s at the beginning of the null 
cone using 



r = r 3 



dr 



1 



dt l 



dr 



2 

N a —1 . 



dt 1 



non+2 



(4.47) 



Notice that this expression requires a value for 8 at level 
n + 2, something that won't be known until the next time 
step. We have found it sufficient to extrapolate from the 
previous time levels using 



n+l 



3/3 



■0 



N 



(4.48) 



Thus, to get v alues for 3 and V we use the matching 
conditions (4.36) and (4.37) along with the above inter- 
polation. For the scalar field, we interpolate the value of 
g from the Cauchy grid, using g — re/). 

The apparent horizon is found on the ingoing null cones 
using the apparent horizon equation which reduces sim- 
ply to V = 0. When the scalar field passes into the 
black hole, the horizon grows outward and we simply 
stop evolving the grid points that are now inside. 



D. Performance 

To evaluate the performance of this approach, we com- 
pare it to a second order accurate, purely Cauchy evolu- 
tion in IEF coordinates, as presented in || . For the com- 
parison shown here, we place the outer boundary of the 
Cauchy evolution at r = 62M and evolve to t = 40M to 
prevent any outer boundary effects from influencing the 
comparison. The scalar field pulse is centered at r = 22M 
and has a mass of 0.5 M. 

FIG H shows the mutual convergence between the 
Cauchy variables in the two codes. This demonstrates 
that the two programs are solving the same problem, and 
provides evidence that the matching approach generates 
the correct spacetime. 
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I 



5 10 15 20 25 30 35 40 
t 

FIG. 6. This plot shows the mutual convergence of the 
Cauchy variables from the matching evolution and the pure 
Cauchy evolution. The evolutions are for a strong, ingoing 
scalar pulse (0.5M) outside a black hole. A convergence factor 
near four indicates second order convergence, while one near 
two indicates first order convergence. 

The reason that the convergence rate appears to drop 
to first order when the scalar field hits the horizon is 
an artifact arising from the motion of the horizon. As 
mass falls into the hole, there is a critical amount which 
causes the horizon to move out by one grid point. If at 
the coarsest resolution the horizon moves a distance dr 
then on the next finer grid it only moves by dr/2, and so 
on. Thus, at different resolutions, the black holes have 
slightly different locations. The resulting shift in the lo- 
cation of the inner boundary causes convergence between 
successively finer numerical solutions to drop from sec- 
ond order to first order. This is an unavoidable diagnostic 
effect due to the comparison of numerical solutions. We 
believe that the numerical solution would converge at sec- 
ond order to an exact solution of the physical problem. 
No exact solutions are known to use for a check of this, 
but for a weak scalar field the horizon does not move and 
we do measure second order convergence throughout the 
evolution. Further, we see the same convergence order 
drop for strong fields in Cauchy-only or null-only evolu- 
tions, and thus are certain it is not due to the matching 
procedure. 



V. CONCLUSIONS 

Our work shows that the matching approach provides 
as good a solution to the black hole excision problem in 
spherical symmetry as previous treatments |^-|]J2?|. It 
also has some advantages over the pure Cauchy approach, 
namely, it is computationally more efficient (fewer vari- 
ables), and is much easier to implement. We achieved a 
stable evolution simply by transforming the outgoing- null 
evolution scheme to work on ingoing null cones, and im- 
plementing it. Achieving stability with a purely Cauchy 
scheme in the region of the apparent horizon is trick- 
ier, involving much trial and error in choosing difference 
schemes. It should be noted, however, that implementing 



the matching may be tricky, especially in higher dimen- 
sions. Whether it is easier than implementing Cauchy 
differencing near the horizon remains to be seen. 

Also, using the ingoing null formulation, we have 
achieved the stable evolution of a Schwarzschild black 
hole in 3-dimensions (the details will be presented else- 
where) and are working on rotating and moving black 
holes. Long term stable evolution of a 3D black hole has 
yet to be demonstrated with a Cauchy evolution. 
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